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Abstract 

Assuming three generations of Majorana neutrinos, we study the different mixing 
patterns that arise from the non-trivial breaking of the flavor group and CP to 
the residual symmetries Z 3 , or Z 2 x Z 2 in the charged lepton and to Z 2 x CP in 
the neutrino sector. All patterns contain only one free parameter 9 and thus mixing 
angles as well as the Dirac and the two Majorana phases are strongly correlated. 
We perform an analytical and a numerical study of all possible mixing patterns. It 
turns out that only four patterns can describe the experimentally measured values 
of the mixing angles for a particular choice of 9 well. All of them predict trivial 
Majorana phases, while the Dirac phase d is maximal for two patterns and trivial 
for the two remaining ones. If 5 is maximal, also the atmospheric mixing angle is 
hxed to be maximal. 



1 Introduction 


By now all three lepton mixing angles have been measured with a certain degree of 
precision (see also iS) 


sin^ ^13 ^ 0 . 0218 ( 9 ) , sin^ O12 ^ 0.304^0:012 > 


23 = 


0.452 

(0.579 


( 1 ) 

for a normal ordering (NO) and in brackets for an inverted ordering (10) of the neutrino 
masses, respectively. Assuming that neutrinos are Majorana particles and none of them 
to be massless, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix Upmns 
that encodes lepton mixing does not only contain the three mixing angles 612 , ^13 and 
023; but also three phases: the Dirac phase 5 and the two Majorana phases a and /3. The 
former can be measured in neutrino oscillation experiments, while a linear combination 
of the latter can be accessible in neutrinoless double beta decay experiments (for a review 
on leptonic CP violation see [^). No direct signals of CP violation have been observed 
in the lepton sector and recent global fits only show a weak indication (below the 3 a 
significance) of a non-trivial value of the Dirac phase 6 [^. 

Many approaches have been pursued in order to describe the data on lepton mixing. 
A particularly promising Ansatz assumes the existence of a flavor symmetry Gf, usually 
finite, non-abelian and discrete, that is broken to different residual groups Ge and G^, 
in the charged lepton and neutrino sectors, respectively (for reviews see [^[^). In (the 
most predictive version of) such a framework all three mixing angles together with 
the Dirac phase 5 can be fixed by the symmetries of the theory Gy, G^ and G^,, if 
the three generations of left-handed (LH) leptons are assigned to an irreducible three- 
dimensional representation 3 of the flavor group G/, i.e. to a representation that cannot 
be decomposed further in Gjj^ We note that this approach does not constrain lepton 
masses and thus all statements made on the predictive power regarding lepton mixing 
angles and the Dirac phase 5 are valid up to possible permutations of rows and columns of 
the PMNS mixing matrix. For neutrinos being Majorana particles surveys of finite, non- 
abelian and discrete subgroups of SU{3) and f7(3), see and references therein, have 
shown that symmetries giving rise to mixing angles which are in good agreement with 
experimental data in general lead to trimaximal (TM) mixing (and thus sin^ 0 i 2 > 
1/3) and a trivial Dirac phase <5 = 0, tt. 

An extension of this approach that involves also CP as symmetry has been proposed 


in (see also 11 as well as [12] ). In this case, also the CP symmetry acts in 
general in a non-trivial way on the flavor space [13] and conditions have to be fulfilled 
in order for the theory to be consistent. The residual groups Gg and Gy are chosen as 
follows: Ge is an abelian subgroup of G/ with three or more elements and Gy is the 
direct product of a Z 2 group contained in Gj and the CP symmetry. Thus, the form of 

^Since lepton mixing like quark mixing only regards LH fields, we do not need to specify in this 
approach how right-handed (RH) charged leptons, and possibly RH neutrinos, transform under G/, 
unless we construct an explicit model in which this approach is realized. 
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these groups is similar to the one in the approach without CP. The advantages of the 
extension with CP are threefold: Majorana phases are also predicted, the Dirac phase 
does not need to be trivial, if the lepton mixing angles are accommodated well, and the 
mixing pattern contains one free parameter 6 . The latter allows for a richer structure of 
patterns that are in good agreement with the experimental data. In particular, mixing 
does not need to be TM. At the same time, being governed by only one free parameter 
all mixing parameters are strongly correlated. Such correlations can be testable and/or 
distinguishable at future facilities 14 . The value(s) of this free parameter that admit(s) 


a reasonable agreement with the experimental data is (are) not fixed by the approach 
itself, but has (have) to be achieved in a concrete model (see for several successful 
models with different symmetries Gj and CP). This approach has already been studied 
for a variety of flavor symmetries: A 4 and ^4 (^[^, A(48) [^, A(96) as well as 

A 5 as flavor group. This group has already been 


A(3n^) and A( 6 n^) with general n 
Here we would like to consider 


employed as flavor symmetry |20| - |23| . In particular, it has been shown to give rise to the 
so-called “golden ratio” (GR) mixing pattern, sin^ @23 = 1/2, @13 = 0 and tan 0 i 2 = !/</ 
with (p = {1 + \/5)/2 Ri 1.618 so that sin^ @12 ~ 0.276|j Very recently, predictions 
of CP phases have been discussed in a scenario with A 5 as flavor group and a CP 
symmetry [26] . Since the authors assume a Klein group and a CP symmetry to be 
preserved in the neutrino sector, the mixing angles are fixed to the values of the GR 
mixing pattern, while possible values of the two Majorana phases depend on the GP 
transformation that is preserved. The latter is not constrained to correspond to an 
automorphism of the flavor group A 5 and thus results obtained in 
In 


26 differ from ours. 


27 


an Ansatz has been pursued in which two GP symmetries are present as residual 
symmetries in the neutrino sector. The combination of these two leads in general to 
a symmetry acting on the flavor space only. Under certain conditions this can be a 
transformation belonging to a finite, non-abelian and discrete group. Therefore, if the 
latter is the alternating group A 5 results of our study can also be achieved using the 
Ansatz with two GP transformations in the neutrino sector. 

In the present paper we analyze the scenario with the flavor symmetry A 5 and a 
CP symmetry comprehensively, since we consider all CP symmetries that correspond to 
involutive ‘class-inverting’ automorphisms of A 5 , all possibilities for Ge, i.e. Ge = Z^, 
Ge = Z^ as well as Ge = Z 2 x Z 2 , and all possible Z 2 subgroups of A 5 as residual 
flavor symmetry in the neutrino sector. All these combinations of symmetries together 
with all possible permutations of the rows and columns of the PMNS mixing matrix are 
subject to an analytical and a numerical study. In particular, we perform a analysis 
using the results of the mixing angles from the global fit [^. As outcome we only find 
four patterns that admit a reasonable agreement with the experimental data, i.e. at 
the 3cr level or better, for a particular choice of the parameter 6 . Two of these four 


^Note that a variant of this approach has been considered for Gf — A(6n^) in which the residual 
symmetry in the neutrino sector is a Klein group contained in the flavor group and a CP symmetry 


19 


^Different versions of the GR mixing pattern are known in the literature that lead to different pre 


dictions for the solar mixing angle in terms of the golden ratio [22[|24[[25] . These are based on different 
flavor symmetries. 
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patterns predict a maximal Dirac phase together with maximal atmospheric mixing, 
while the other two ones lead to a trivial Dirac phase and in general non-maximal 023- 
Majorana phases are trivial for all four patterns. Thus, two out of these four lead 
to no CP violation. This fact can be traced back to the existence of an accidental CP 
symmetry, common to the charged lepton and neutrino sectors, as we discuss. As regards 
the reactor and the solar mixing angles, we note that ^13 is in general accommodated 
well, whereas O 12 is subject to non-trivial constraints in all four cases: two of the four 
patterns give rise to a lower bound sin^ 9 i 2 > 0.276, the value of the GR mixing pattern, 
one incorporates TM mixing and thus sin^ 612 ^1/3 and the remaining one entails an 
upper limit sin^ 612 ^ l — ~ 0.345, a value that is associated with a different version 

of the GR mixing pattern [24] . 

The paper is organized as follows: in section we recapitulate the approach with a 
flavor and a CP symmetry as well as the main features of the group A 5 . We also discuss 
the admitted CP transformations and relegate further details regarding their relation 
to the automorphisms of A 5 and their nature to appendix Section contains the 
analytical study of all patterns that can lead to a good agreement with the experimental 
data as well as the results of our analysis. We summarize our main results in section 
1^ Besides appendix we include appendix that contains our definitions of mixing 
angles, CP phases and corresponding CP invariants Jcp, Ii and l 2 - 


2 Approach 


We briefly recapitulate the essential ingredients of the approach and summarize the 
necessary information on the group A 5 . We list the candidates of generators of residual 
flavor symmetries in the charged lepton and neutrino sectors as well as the CP trans¬ 
formations. At the end of this section we also comment on the possible presence of 
an accidental CP symmetry common to charged leptons and neutrinos. Since we focus 
on the case in which the generations of LH leptons are assigned to an irreducible three- 
dimensional representation 3 of the flavor group, the CP transformation X as well as the 
elements of the flavor symmetry are represented by (unitary, complex) three-by-three 
matrices in the following. 

Let us consider a theory with a flavor group Gf = combined with a CP symmetry 
that in general also acts non-trivially on the flavor space 13 . The CP transformation 
A is a unitary and symmetric matrix 


AA^ = AA* = 1 


(2) 


The requirement that A, the CP transformation associated with the CP symmetry pre¬ 
served in the neutrino sector and subject to the condition in Q, must be a symmetric 
matrix has been shown in to be necessary, since otherwise the neutrino mass spectrum 
would be partially degenerate and, consequently, inconsistent with experimental obser¬ 
vations . In order to ensure a consistent combination of the flavor and CP symmetry 
we require that the subsequent action of the CP transformation, an element of the flavor 
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group and the CP transformation is equivalent to the action of an (in general different) 
element of the flavor group 


{X-^AXf = A', 


( 3 ) 


10 


12 


a CP 


with A and A' representing (different) elements of As shown in 

transformation corresponds to an automorphism of the flavor group. In particular, our 
request that X fulfills Q renders this automorphism involutive. Following the discussion 
in 11 this automorphism should be class-inverting, i.e. the image of the element g under 


the automorphism has to lie in the same class as the inverse of g. This is guaranteed 
for the three-dimensional representation 3 by the fulfillment of the condition in Q. As 
we show in appendix [A| the automorphisms corresponding to the CP transformations we 
consider in our analysis are also class-inverting when acting on the other representations 
of the flavor group A^. 

The residual symmetry in the neutrino sector is assumed to be the direct product of 
a Z 2 symmetry contained in the flavor group and the CP symmetry. Thus, the matrix Z 
representing the generator of the former symmetry and the CP transformation X have 
to fulfill 


XZ* - ZX = 0, 


( 4 ) 


which is a particular case of the condition in ([^ . We note that the presence of the residual 
symmetries given by Z and X implies the existence of a second CP transformation 
Y = ZX in the neutrino sector that fulfills the same conditions in ([^-Q as the CP 
transformation X [^. 

In the charged lepton sector, in contrast, we take as residual group an abelian sub¬ 
group of the flavor symmetry that offers the possibility to distinguish among the three 
generations of charged leptons, i.e. this group has to have at least three different ele¬ 
ments. The residual group can be described with a set of generators Qi, z = 1, 2,... that 
commute. 

The derivation of lepton mixing in this scenario has already been presented in detail 
in and we only mention it briefly. In the charged lepton sector, the residual group 
generated by Qi constrains the charged lepton mass matrix mi , here given in the right-left 
basis, 

Q\ m\mi Qi = m\mi . (5) 

The matrices Qi are diagonalizable by the unitary matrix C/g, i.e. ulQiUe is diagonal, 
which is determined up to the ordering of its columns and possible overall phases of the 
single columns. As a consequence, also the combination 


Ul m\mi Ue is diagonal . 


( 6 ) 


Thus, Ue diagonalizes the charged lepton mass matrix as regards LH charged leptons. 
In the neutrino sector, the light neutrino mass matrix miy is subject to the following 

^We use throughout this paper lowercase letters for the abstract elements of the flavor group A 5 and 
capital letters for the matrix representatives (in the representation 3). 
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condition^ 


niiy Z = niiy and X niu X = m* , 


( 7 ) 


if the Z 2 and CP symmetry are imposed. Without loss of generality we can choose a 
basis such that 


x = nn'^ and zn = dmg{{-iy\i-iy\{-iy^), (8) 

with n being unitary and Zi = 0,1. Since Z generates a Z 2 symmetry, two of the three 
parameters Zi have to coincide. The matrix combination ff" my Q is then real and block- 
diagonal and thus can be diagonalized by a rotation Rij{9) in the (ij)-plane through an 
angle 9 that is determined by the matrix entries of my Q. Their actual values are 
in general not predicted in this approach and thus 9 is taken to be a free parameter in 
the interval between 0 and tt in the following. The (ij)-plane is fixed by the degenerate 
sub-sector of 01" Z Q, i.e. the two Zi and Zj that are equal. The positiveness of the 
light neutrino masses (a vanishing neutrino mass can also be included) is ensured by 
the diagonal matrix Ky with entries ±1 and ±i on its diagonal. So, the contribution to 
lepton mixing from the neutrino sector is given by 


Uy=nRiy9)Ky, 


( 9 ) 


and the PMNS mixing matrix resulting from this approach reads 

UpMNS = U^e^Rij{9)Ky. 


( 10 ) 


It is important to note that this mixing matrix is only determined up to permutations 
of its rows and columns (and unphysical phases), since this approach does not make 
any predictions concerning the mass spectrum of charged leptons and neutrinos. For 
example, if not embedded in a model context, see e.g. first reference in 15 , one cannot 
predict whether neutrinos follow NO or 10. It is also worth to emphasize that a mixing 


matrix of the form in ( 10 ) has one column that is determined by group theory only and 


that does not depend on the free parameter 9. 

We note also that two tuples {Q,Z,X) and {Q',Z',X') lead to the same physical 
results, if the generators of the symmetries are related by a similarity transformation 


n^Q'n = Q , z'n = z and x'n* = x. 


( 11 ) 


The alternating group A 5 describes the even permutations of five distinct objects. It 
is isomorphic to the icosahedral rotation group I. It has 60 different elements, organized 
in five conjugacy classes, and, thus, possesses five irreducible representations: 1, 3, 3', 
4 and 5. All representations apart from the singlet are faithful]^ The group A 5 can be 

®We do not need to specify the generation mechanism of neutrino masses unless we construct an 
explicit model. So, this mass matrix can arise from integrating out heavy RH neutrinos, from Higgs 
SU{2)l triplets acquiring a vacuum expectation value, etc.. 

® A representation is called faithful, if all elements of the group are represented by different matrices in 
this representation. So, in the case of A 5 the 60 group elements are represented by 60 different matrices 
in the representations 3, 3', 4 and 5. 


5 





generated with two generators s and t that fulhll the relations 

= e , = e and (st)^ = e , 


( 12 ) 


with e denoting the neutral element of A^. Since we would like to assign LH leptons to 
triplets, we are particularly interested in the representations 3 and 3k The explicit form 
of the generators s and t in the representation 3 can be chosen as 



and T = 


0 

A $ 


0 

0 


with 


= l{l + V5) 


and 


a. 27r 
= — . 
5 


(13) 


(14) 


The generators in 3' are easily obtained from S and T in (13) by using the combination 
T'^ ST^ S and as generators, see also 22 . This shows immediately that the set of 


all matrices describing the representations 3 and 3' is the same and thus all conclusions 
obtained in a comprehensive study of mixing using the representation 3 also hold for 3k 
Consequently, it is irrelevant for our analysis whether LH leptons are in 3 or 3^ of A 5 
and, without loss of generality, we assume in the following that LH leptons transform as 
3 of A5. 

The group A^ has several subgroups. In particular, the group contains 15 elements 
that generate a Z 2 symmetry which give rise to five distinct Klein groups. These are, 
like the Z 2 generating elements themselves, all conjugate to each other. Since we make 
explicit use of these elements we mention them here 


f 1 = S , f2 = st'^st^st'^ , ^3 = , ^4 = t'^st , f5 = st^st^s , 

stasis , Vj = tst‘^ , V% = st^st^s , Vq = stst^st^ , IJlo = st^st , 

vii = t^st^ , vi 2 = tst^st^s , ni3 = tst^s , rJi4 = t^st^ , vi^ = st^st^st . 

They form the Klein groups Ki 

Ki = {vi,V 2 , f3, e} , K 2 = {r’4, t’5, '^6, e} 

K3 = {U7, Us, Ug, e} , Ki = {uio, Uii, U12, e} and = {uis, U14, U15, e} , (15) 

see also [22] . Furthermore, there are ten Z^ and six Z5 subgroups. Also these are all 
conjugate to each other. For the complete list of generating elements of these subgroups 
see appendix B in [^. So, three different types of groups, Z3, Z5 and Z 2 x Z2, can 
function as residual symmetry Ge in the charged lepton sector. 

The form of the CP transformations we consider is 


X = I/Ao (16) 

with Xq 

1 0 0 \ 

001 (17) 

0 10 / 
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and V is the matrix representative of a Z 2 generating element, i.e. = 1, or V = 1. So, 
we find in total 16 different possible CP transformations. All of them fulhll Q. Notice 
that the CP transformation Xq gives rise to the so-called fi-r reflection symmetry 
whose phenomenological consequences have been studied in detail in the literature. 

As one can check, it holds for Xq 


28 


(Ap-^SAo)* = S and (Ag-^TAo)* = T . (18) 

Thus, for the generators S and T in the representation 3 the condition in Q is valid for 
A = A'. As a consequence, the corresponding automorphism is the trivial one 


s ^ s and t —>• t . 


(19) 


The CP transformations X = V Xq correspond to different inner automorphisms of 
the group A5 (i.e. automorphisms whose action on the elements of the group can be 
represented by a similarity transformation with a(nother) group element) that map the 
generators s and t in the following way 

s —)■ vsv~^ and t —> vtv~^ . (20) 

The automorphism group of A5 is the symmetric group S'5 and the group of inner au¬ 
tomorphisms is isomorphic to A5 itself. As we show in appendix [A| the 16 different CP 
transformations we consider correspond to the 16 class-inverting involutive automor¬ 
phisms of A5. We thus discuss all CP transformations that can be consistently imposed 
according to |llj and that fulfill the requirement in Q. 

A last condition that needs to be examined is the constraint in (Q, namely whether 
the Z 2 generator commutes with the chosen CP transformation. We hnd that for each 
Z four different possible CP transformations A are admitted: for Z being one of the 
non-trivial elements of the Klein group Aj, i = 1,..., 5, the CP transformation X = V Xq 
with V belonging to the same Klein group Aj (this time K = 1 is included) is a viable 
choice. Taking into account that the Z 2 generator Z and the CP transformation A 
automatically imply the existence of a further CP transformation Y, Y = ZX, we can 
reduce the number of independent choices of CP transformations for each Z 2 generator 
Z to two. 

Eventually, we mention that it can happen that an accidental CP symmetry Y, 
common to the charged lepton and the neutrino sectors, exists that leads to trivial 
CP phases. This can be checked by searching for a transformation Y that fulhlls the 
following constraints 

Y* mjrniY = {mjmiy and Y m,yY = m* . (21) 

These are equivalent to the conditions involving the generators of the different symme¬ 
tries 


Q^Y-YQf = 0, (22) 

ZY - YZ* = 0 and XY* - YX* = 0 , (23) 

with Y being a diagonal and real matrix in the neutrino mass basis, i.e. Y U* is 
diagonal and real. For details see [^. 
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3 Lepton mixing 


In order to study mixing comprehensively, we analyze all possible combinations of resid¬ 
ual symmetries in the charged lepton and neutrino sectors, i.e. all possible Gg and Gi^ 
{Z 2 generators and CP transformations X). These can be expressed as tuples of gen¬ 
erators Qi of Gg and Z and X for Gy. As mentioned in section Gg can be either 
a Z3, Z5 or a Klein group, while one of the 15 different Z 2 symmetries that each can 
be consistently combined with four different CP transformations X specify the residual 
group Gy. Instead of computing the mixing pattern for all these we can greatly reduce 
the number of cases that we need to study by applying similarity transformations as 
well as the fact that for a pair Z and X also the pair Z and Y = ZX leads to the 
same mixing pattern. We allow for all possible permutations of rows and columns of 
the mixing matrix. Furthermore, we consider both possible neutrino mass orderings NO 
or 10 in our (numerical) analysis. We exclude all patterns that cannot accommodate 
the experimental data on lepton mixing angles at the 3 a level or better for (a) certain 
choice(s) of the free parameter 9. As a consequence, we end up with in total only four 
cases. These we call in the following Case I to III and Case IV-PI/Case IV-P2. We 
first study the mixing patterns analytically for each possible residual symmetry Gg in 
the charged lepton sector and then show the results of a analysis of these patterns, 
since they can accommodate the experimental data best. We briefly comment on a fifth 
mixing pattern that can fit the data also well apart from the solar mixing angle whose 
value turns out to be slightly smaller than the lower 3 a bound on sin^ 612 . 

3.1 Gg = Z 5 : Case I and Case II 

If we consider as Gg a Z5 symmetry, we find six different categories of tuples (Q, Z, X) 
with Q being a generator of a Z5 group, taking into account the mentioned operations in 
order to relate different tuples {Q,Z,X). We can show that for each of these categories 
we can hnd a representative tuple {Q,Z,X) with Q = T. Thus, the mixing matrix Ug 
resulting from the charged lepton sector is the unit matrix, up to permutations of its 
columns and unphysical phases. Out of these six representative tuples four lead to a 
mixing pattern that we dismiss, because the column that is fixed by group theory is not 
compatible with the data at the 3 a level or better [^. The two remaining representatives, 
which we can choose as 


(Q, Z, X) = (T, T^ST^ST^, SXq) (Case I) (24) 

and 

(Q, Z, X) = (T, ST^ST, Xo) (Case II) , (25) 

give rise to a mixing matrix with a column whose components have the absolute values 

/ sin If \ / 0.526 \ 

cosf/V 2 ^ 0.602 , (26) 

\ cos(/?/-v/2 / \ 0.602 / 
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where we defined for convenience 


iw.(p = !/(/). 


(27) 


This column can only be identified with the second one of the PMNS mixing matrix, if 
good agreement with the experimental data should be achieved. Notice that the ordering 
of the components in this column as well as its position within the PMNS mixing matrix 
are not fixed by the approach we are using, since no constraints on the lepton masses 
are imposed. 

We first derive the mixing pattern for the tuple in (24). As explained Ue = 1 and we 
can take rij to be 



—\/2i siny? 
i cos ip 
i cos ip 



(28) 


that fulfills ([^ for X and Z chosen as in (24). Since zi and of the (diagonal) com¬ 
bination nj Zrij are equal, see (8) for definition, the correct indices ij of the rotation 
matrix Rij{9) in @1 are ij = 13. Thus, the PMNS mixing matrix reads 


UpMNS = -^ 13 ( 6 *) Ku 


(29) 


We can extract the mixing angles from (29) in the usual way and find 

2 . o „ 1 

10 


sin^ 6*12 = 


2 -|- (3 -|- \/5) cos^ 9 


sin^ 013 = 


5 -|- Vo I sin^ 9 


1 


sin"' 023 = X - 


y^2(5-h\/5) sin 20 
7 -k \/5 -k (3 -k \/5) cos 20 


(Case I) 


(30) 


Furthermore, we can derive the following exact sum rule among the solar and the reactor 
mixing angles 

2 „ sin^ ip 0.276 


sin 012 = 


> 0.276 . 


(31) 


1 — sin^ 013 1 — sin^ 0i3 

This sum rule can also be directly obtained from |f7e2p = sin^ (/j = 1/(1 -|- cj)'^) ~ 0.276. 
Using for sin^ 0i3 its best fit value (sin^ 0i3)^^ = 0.0219 we find for the solar mixing angle 
sin^ 012 ~ 0.283 which is within the 3cr range, see ([^. This value coincides with the one 
obtained from a analysis, see table[^ The non-trivial lower bound on the solar mixing 
angle, see ( |3l| ), is also nicely seen in figure (and implicitly also in the sin^ 0i3-sin^ 0i2 
plane in figure . For the atmospheric mixing angle a simple approximate relation to 
the reactor mixing angle is given by 


sin^ 023 ~ - ^1 ± (1 — Vs) sin 013^ « 0.5 =F 0.618 sin0i3 ss | 


0.409 for 0 < 7r/2 
0.591 for 0 > 7r/2 


(32) 


if we use again the best fit value (sin^ 0i3)'^^ = 0.0219. The first subleading term arises at 
sin^ 013 and can thus be safely neglected in this estimate. We can clearly see having a look 
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at the different symbols in figures and that values of 0 < 7r/2 lead to sin^ 023 <1/2, 
while larger values, 7r/2 < 0 < tt, entail sin^ 023 > 1/2. This is also confirmed by the 
two different ‘best fitting’ values 0bf 0.174 and 0bf ss 2.967 that are obtained for NO 
and 10, respectively, see table 

Note that the formulae for the mixing angles 0i2 and 0i3 remain invariant, if we 
replace 0 by vr — 0, while the relative sign in the expression for sin^ 023 changes. The 


same effect can be achieved, if we consider the PMNS mixing matrix in (29) with second 
and third rows exchanged. 

As one can check, all CP invariants Jcp, Ii and I 2 vanish exactly and thus an 
accidental CP symmetry must be present. Indeed, there is one, namely 



(33) 


that fulhlls the conditions in (21)-(23) for the tuple {Q,Z,X) shown in (24). The van¬ 
ishing of sin (5 can also be confirmed by verifying that the condition presented in [27] is 
fulfilled. 

In a similar manner we can study the lepton mixing that can be obtained for the 
second tuple (Q, 2', X), the one in (25). The form of the matrix IIjj that fulfills ([^ for 
Z = ST'^ST and X = Xq can be chosen as 


—\/2cos(/? — \/2 sin(^ 


fill = — I - 


3— 3i$ 


v /2 


sm(f 


cos (p 


_ 2 7 • _ 2 7 <J> 

—e simp e cos ip 


0 

-e^ 
5-3i<I>/4 


(34) 


Like in the preceding case, also here the necessary rotation Rij{0) is in the (13)-plane. 
Thus, taking into account that Ue is trivial, the PMNS mixing matrix is of the form (up 
to permutations of rows and columns and unphysical phases) 


UpMNS = RisiO) Ky 


(35) 


The predictions for the solar and the reactor mixing angles are the same as for the first 
tuple, namely 


sin^ 012 = 


2 -h (3 -h \/5) cos2 0 


sm 


9 i3 = — ( 5 -|- V5 ) sin^ 0 (Case II) . 


(36) 


Thus, also the sum rule and the estimate given in and below (31) hold. In figure 
the non-trivial lower bound on sin^ 0i2 is visible also for Case II. The results for the 
atmospheric mixing angle and the Dirac phase instead are different, since both of them 
are predicted to be maximal 


1 


sin^ 023 = - and | sindl = 1 (Case II) 


(37) 
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The sign of sin <5 depends on whether 9 ^ 7 r/ 2 . The corresponding Jarlskog invariant 
Jcp read^ 


Jcp 


1 

20\/2 


\lb + Vh sin 29 . 


(38) 


Like in the first case also in this case the Majorana phases are trivial. We can check 
that the conditions found in |27| for cos J = 0 and 023 = 7 r /4 are fulfilled in the case at 
hand. Furthermore, we can verify with the help of the formulae given in that both 
Majorana phases must be trivial. 

We notice that the replacement of the parameter 0 by vr — 0 does not change the form 
of the mixing angles, while the sign of Jqp and, consequently, of sin <5 changes. The very 
same result is achieved, if we exchange the second and third rows of the mixing matrix 
in (35). Consequently, we expect to find two best fitting values of the parameter 0. This 
is confirmed by the results of the analysis in table As expected, the sum of these 
two best fitting values equals vr. 

One of the categories initially dismissed by the criterion that the absolute values of 
the group theoretically fixed column of the PMNS mixing matrix should agree at the 
3(T level or better with the values given in might still be interesting in a concrete 
model in which (small) corrections can lead to the agreement with experimental data. A 
representative of this case is the tuple (Q, Z, X) = (T, S', Xq). The absolute values of the 
elements of the fixed column are (cos<y 9 , sin(/?/\/ 2 , sin ip/V^) PS (0.851,0.372,0.372)^. 
Thus, this column can be identified with the first one of the PMNS mixing matrix. This 
pattern fails to describe the data well without corrections mainly because of the tight 
relation between the solar and the reactor mixing angle that can be derived. We find 
sin^ 012 = 1 — 10 ~ 0.276 that leads for (sin^ 013)'’^ = 0.0219 to a too small 

solar mixing angle sin^ 0 i 2 ~ 0.260. At the same time, the atmospheric mixing angle is 
maximal. The Dirac phase is also maximal, whereas both Majorana phases are trivial. 
So, this case shows strong similarities to Case II with the representative tuple shown in 


(25). 


3.2 Ge = Zs’. Case III 


If we do the same analysis for the case Ge = Z 3 , we find eight categories of tuples 
{Q,Z,X). We can always fix Q = T'^ST'^. In this case Ue is not trivial anymore and is 
of the form 


Ue = 


/ 

V 


7+3\/5 

3(5+V5) 

2 


V3(5+C5) 

2 


V 15 

- 3 ^ \/75 + 30v/5 
i ^ ^75 + 30\/5 


2 ^ 

V3(5 +a/5) 

-^ + ^\/75 + 30v/5 


(39) 


^For sin 26 = 0 the Jarlskog invariant vanishes. If this happens, one of the mixing angles becomes 
either 0 or 7r/2 and thus the Dirac phase 5 becomes unphysical. Clearly, these values of 6 are highly 
disfavored by experimental data. 
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Using the representatives of the eight different categories we see that indeed only one of 
these can lead to a mixing that is compatible with experimental data. The column that 
is fixed by group theory is in this case TM 


1 

71 



0.577 \ 
0.577 
0.577 / 


(40) 


and has therefore to be identified with the second column of the PMNS mixing matrix. 
Immediately, we know that the solar mixing angle has a lower bound, sin^ 0i2 ^1/3, 
see the sum rule in (46), the result of the analysis in table and the lower bound in 
figure We use as representative the tuple 


(Q, Z, X) = {T^ST'^, ST^ST^S, Xq) (Case III) 
A possible admitted form of the matrix U is 


(41) 




\/2cos(/? 0 —\/2sin<^ 


7 ^ 

e sm(f 


_ 7 ^ 

e siiKf e 


e" ^ cos if 


(42) 


cos ip 


As can be seen, the form of the matrix is quite similar to the ones used in the other 
cases, see (28) and (34). The matrix Uu is composed as follows 

u, = niiiRi3{e)K,, (43) 

since, as in the cases above, zi and Z 3 of the matrix combination Ujjj-ZQjp are equal. 
The PMNS mixing matrix is then given by 


UpMNS = ul Hip i?l3(6') Kj, . 


(44) 


We can extract the following results for the solar and the reactor mixing angles 


sin^ 012 = 


„ , and sin^ 0i3 = ^ (1 - sin20) (Case III) 

2 + sm 20 4 


that fulhll the exact - and well-known - sum rule 

sin^ 012 = 


3(1 - shr 013) 


>1 


(45) 


(46) 


If we use (sin^ 013 )'’^ = 0.0219 as best fit value for sin^ 0i3, we arrive at sin^ 0i2 ~ 0.341 
which is the value that is also obtained in our x^ analysis, see table The atmospheric 
mixing angle as well as the Dirac phase are, as in Case II, maximal 


sin^ 0' 


23 = - and I sin 51 = 1 (Case III) 


(47) 
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Again, the actual sign of sin <5 depends on the choice of 6 . 
invariant is 


cos 26 


The form of the Jarlskog 


(48) 


If cos 20 vanishes, Jqp equals 0. Since at the same time one mixing angle becomes 0 or 
7 r/2, the Dirac phase 5 turns out to be unphysical for cos 20 = 0. As happened for Case 
I and II, also here the mixing pattern that is well compatible with experimental data 
gives rise to trivial Majorana phases. 

We note that the formulae of the mixing angles remain invariant, if we replace 0 with 
7 r/2 — 0. Thus, we expect (at least) two best fit solutions. This expectation is confirmed 
by our analysis, see table Replacing 0 with 7r/2 — 0 leads, at the same time, to an 
additional sign for Jcp and thus sin (5. Similarly, the exchange of the second and third 
rows of the PMNS mixing matrix in (44) does not alter the results for the mixing angles, 
but changes the sign of the Jarlskog invariant and thus of sinJ. 


3.3 Gp = Zo X Zn'. Case IV-Pl and Case IV-P2 


For the remaining possibility = Z 2 x Z 2 we find that all admitted combinations of 
Qi, Q 2 with Z and X describing the residual symmetry Gy = Z 2 X GP can be classified 
in four different categories of tuples {{Qi,Q 2 } , Z, X). Thus, it is sufficient to calculate 
the mixing pattern for one representative of each category. Note that we can always 
choose a representative for which Gg = Ki, see (15), i.e. Qi and Q 2 can be chosen as 
Qi = S and Q 2 = T'^ST^ST'^. So, the form of the matrix Up is 





(49) 


It turns out that only one category of tuples is capable of accommodating the experi¬ 
mental values of the mixing angles well for a particular choice of the parameter 0, while 
the other three ones fail to do so. In particular, two out of these three lead to patterns 
with only one non-vanishing mixing angle, since the generators Qi are diagonalized by 
the same matrix as the Z 2 generator Z. A representative of the category that allows for 
good agreement with the data is 

{{Qi,Q 2 },Z,X) = {{S,T‘^ST‘^ST‘^],ST‘^ST,Xq) (CaselV). (50) 

We can check that the column that does not depend on the free parameter 0 has com¬ 
ponents with absolute values of the form 


/ </- \ / 0.809 \ 

- 1 U 0.5 . (51) 

\l/ct> ) \ 0.309 ) 

Thus, this column can be identified with the first one of the PMNS mixing matrix. We 
call this situation Case IV-Pl. We note that we can exchange the second and third 


13 





components of the vector in (51), i.e. we can exchange the second and third rows of the 
resulting PMNS mixing matrix, and also obtain good agreement with the experimental 
data. This situation is denoted by Case IV-P2 in the following. The crucial change 
occurs in the predicted value of the atmospheric mixing angle, see (56) and (60). 

for Z = S'l^T and X = X 


A unitary matrix that fulfills the conditions in 
siny? —\/2 


0 IS 


1 


COS^ ip 


V 2 


cos ip sm p> 


Q.IY =I e *^/^cos (/9 e + cos y?) sin(/? e cos — sin^ (/?) 


_ 2 ? 

—e cos(/9 


-2i4> ^pl5i$/4 _ 


COS 




sm(^ 


e cos (/9 + sin^ (/?) 


Since the diagonal matrix Hjy Z Hjv reveals equal Z 2 and Z 3 in the convention of nSn , 
the form of the neutrino mixing matrix is given by, up to permutations of columns and 
unphysical phases, 

= Qiy R23{e) . (53) 

Taking into account the contribution to leptonic mixing coming from the charged lepton 
sector that is encoded in the matrix Ue in (49), the PMNS mixing matrix is of the form, 

(54) 


up to permutations of rows and columns and unphysical phases, 

UpMNS = ^IV R2'i{G) Kv ■ 

We hnd for the mixing angles the following expressions 


sin^ 012 = 


sin^ 023 = 


(5 — 2 1 / 5 ) cos^ 0 

1 + (5 - 2 v/5) cos2 0 ' "" 8 

3 (5 - \/5) + (9 - 5 \/5) cos20 + 8sin20 


sin^ 013 = X (5 — 'v/s') sin^ 0 (Case IV) (55) 
8 V / 


(Case IV-Pl) 


(56) 


25 - 3 + 5 (3 - v/5) cos 20 

For this case, as mentioned above, the absolute values of the elements of the first column 
of the PMNS mixing matrix are ordered in the same way as in (51). From (55) we can 
derive an exact sum rule relating the solar mixing angle to the reactor one 

3 + \/5 


sin 012 = 1 — 


0.331 


(57) 


8(1 -sin^ 013 ) 

if we insert the experimental best fit value (sin^ 013)*^^ = 0.0219. We note that in this case 
a non-trivial upper bound on sin^ 0 i 2 exists, namely sin^ 0 i 2 < 1 - (p^/4 PS 0.345. This 
can also be directly derived from the constraint that |?7eip = 0^/4 = (3 + \/5)/8 ~ 0.655. 
This bound is marked with a (dashed red) vertical line in hgure[^ Furthermore we can 
obtain an approximate relation of the atmospheric mixing angle and sin 0 i 3 


sin^ 0 ' 


23 


1 

10 


5 - V5 + 


Y^5 + 2\/5 sin013 « 0.276 + 1.23 sin0i3 , 


(58) 


using that 0 lies in the interval [0,7r/2] (and thus cos0 = -\/l — sin^ 0)lu Subleading 


(59) 


corrections are of order sim 0i3 at maximum. For (sin^0i3)“^ = 0.0219, we find 

sin^ 023 PS 0.459 . 


® For 7r/2 < 6 < tt we can derive a similar relation that shows, however, that the measured values of 
the atmospheric and the reactor mixing angle cannot be accommodated well at the same time. 
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This estimate is consistent with the result of our analysis, see table 

If we permute the second and the third rows, the form of the reactor as well as of the 
solar mixing angle is the same, while the atmospheric one turns out to be 1 — sin^ 023 ) 
i.e. here the atmospheric mixing angle reads 


sin^ 023 = 1 - 


3 (5 - \/5) + (9 - 5 ^/5) cos20 + 8sin20 
25 - 3 \/5 + 5 (3 - \/5) cos 20 


(Case IV-P2) . 


So, in this case the approximate sum rule relating 623 and 613 is given by 


(60) 


sin2 023 ~ — f5 + \/5) - - V5 + 2\/5 sin0i3 « 0.724- 1.23 sin0i3 . (61) 

10 V / 5 

Again, we assume 0 < 0 < 7r/2 (and for 0 lying in the interval [7r/2,7r] see footnote]^. 
For (sin^ 013)*’^ = 0.0219 we get 

sin^ 023 0.541 (62) 

which is consistent with the value obtained from the fit) see table 

All CP phases are trivial. This points towards an accidental CP symmetry in the 
theory. This is clear, since the CP transformation X = JTq is not only present in the 
neutrino sector, but also - as one can check explicitly - in the charged lepton one. 

Note that there is no evident symmetry as regards the parameter 0 in the formulae 
of the mixing angles: while 0 i 2 and 0 i 3 remain invariant, if we replace the parameter 0 
with TT — 0, this is not the case for the atmospheric mixing angle and thus we expect 
in general only one value of 0 for which the mixing angles can be accommodated best. 
This is confirmed in our numerical analysis, see table 


3.4 Numerical discussion 


In the following we present our results of a analysis for the different cases. Case I 
through IV-P2. The x^ function is defined in the usual way 


with 


2 2,2,2 

X — X12 + X13 + X23 

2 f sin^ 6 ij — (sin 


Xij 


(sin^ Ojj) 


,bfN 2 




for ij = 12,13, 23 


(63) 

(64) 


sin^ 9ij are the mixing angles derived in the different cases, e.g. see (30) for Case I, that 


are the best 


depend on the continuous parameter 0, ranging from 0 to vr, (sin^ Oij) 
fit values and Ujj the la errors reported in (Q. Note that these errors also depend on 
whether sin^ 9ij is larger or smaller than the best fit value. Since the global fit results 
for the mixing angles (slightly) differ for the case of NO or 10, we consider these two 
separately and calculate for all patterns the x^ function Xno under the assumption of 
NO being realized in nature and Xio I®. In particular, in doing so we do not take 
into account the fact that NO is slightly disfavored by Ax^ = 0.97 compared to 10 [^. 
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Case I 

Case II 

Case III 

Case IV-Pl 

Case IV-P2 

{Qi,z,x) 

{T,T^ST'-^ST^,SXo) 

{T,ST^ST,Xo) 

(T^ST^,ST^ST^S,Xo) 

({S, T^ST^ST'^] ,ST'^ST,Xo) 


NO 

10 

NO 

10 

NO 

10 

NO 

10 

NO 

10 

Amin 

5.64 

3.46 

4.04 

7.74 

8.84 

12.56 

4.48 

11.80 

6.19 

6.43 

^bf 

0.174 

2.967 

f 0.175 
\ 2.967 

f 0.604 
\ 0.967 

f 0.603 
\ 0.967 

0.254 

0.258 

0.255 

0.254 

sin^ 012 

0.283 

0.283 

0.283 

0.283 

0.341 

0.341 

0.331 

0.330 

0.331 

0.331 

sin^ 013 

0.0217 

0.0219 

0.0218 

0.0220 

0.0217 

0.0218 

0.0219 

0.0225 

0.0220 

0.0218 

sin^ 023 

0.408 

0.592 

0.5 

0.5 

0.475 

0.478 

0.524 

0.525 

Jcp 

0 

0 

=F0.0325 

=F0.0326 

±0.0342 

±0.0342 

0 

0 

0 

0 

sin (5 

0 

0 

Tl 

±1 

0 

0 

0 

0 


Table 1: Results of analysis for Case I through Case IV-P2, displayed separately for the 
assumption of NO and 10. We remind the reader of the representative tuples {Qi,Z,X) that we 
have chosen for the different cases. smallest value of that can be obtained for 

a particular mixing pattern at the best fitting value (s) We mention for each pattern and 

neutrino mass ordering only the lowest value of that can be achieved. All values of sir? 9ij 
are obtained at the given O^j^. For Case II and Case III that predict maximal atmospheric mixing 
and a maximal Dirac phase S, the sign of the Jarlskog invariant Jcp (and of Sind) depends on 
the chosen value of and upper (lower) signs correspond to the smaller (larger) value of Sjy. 
Majorana phases are trivial in all cases, sin a = 0 and sin/3 = 0. 


A mixing pattern is considered to agree reasonably well with the experimental data, if 
< 27 and all mixing angles sin^ 0*^ are within the 3 a intervals given in The 

functions x^q and Xjq are minimized at the best fitting point (s) 9 = 0bf and we 
only report the global minimum/a in table for each case Since the indication of a 
preferred value of the Dirac phase 6 coming from global fit analyses is rather weak [^, 
i.e. below the 3 a significance, we do not include any information on 6 in the function 
in (63). 

Our findings for the different cases are summarized in table As one can see, these 
results agree well with our analytical estimates and observations made in subsections 
3.1 3.3 In particular, the sum of the two best fitting values of 9 for NO and 10 in Case 


I) ^bf, NO ^bf, lO) approximately equals vr, since the formulae for the mixing angles 
012 and 013 are invariant under the transformation 0 —)• vr — 0 , while sin^ 023 turns into 
cos^ 023. Similarly, the sum of the two best fitting points 01,11 and 0 i,f _2 ((almost) the 
same for NO and 10) equals tt in Case II. Also related to the symmetry properties of 
the formulae for the solar and the reactor mixing angles is the observation in Case III 


® The upper limit on X^i X^ 27, is chosen, since it results from summing three 3 ct gaussian errors 
for one degree of freedom. 

We have tested the validity of our x^ analysis by constructing a likelihood function to fit the various 
cases that uses the one-dimensional projections provided in . These results are consistent with those 
in table up to the fact that the roles of the local and the global minimum in Case I become exchanged 
for NO. However, the difference between these two minima turns out to be statistically insignificant. 
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Figure 1: Results for the mixing angles sin^ 023 sin^ 0i2 for Case I to III and Case IV-Pl 
and IV-P2. The plot on the left side assumes NO, while the one on the right 10. The curves 
of the various cases are superimposed to the experimentally preferred 1 cr (green), 2 a (yellow) 
and 3(7 (blue) areas at two degrees of freedom adapted from 0 In drawing these areas we have 
subtracted Ax^ = 0.97 in the case of NO so that the minimal value of Ax^ is zero for both mass 
orderings. The experimental best fit value is indicated with a cross. The best fitting value of 6 
is indicated on each curve with a red dot. For clarity we mark on the curve belonging to Case 
I some particular values of the parameter 9 with different symbols in black: star for 9 = tt/ 19, 
square for 9 = T:fl7, dot for 9 = IGtt/I? and triangle for 9 = IStt/IQ. In addition, an arrow 
indicates on each curve the direction of increasing 9 in the interval 0 < 9 < 9yj-^ (the arrow 
always belongs to the closest label I, II, etc.). Note that the curves belonging to Case II and Case 
III partly overlap. Dashed red vertical lines indicate the non-trivial lower (upper) bound on the 
solar mixing angle in Case I, II and III (Case IV-Pl and IV-P2). 


that the two best fitting points (for NO and 10) and 25 sum up to 7 r/ 2 . Case IV 
does not reveal such a symmetry in the parameter 9 and thus we discuss in this case the 
results corresponding to two different permutations, Case IV-Pl and Case IV-P2, that 
are related by the exchange of the second and third rows of the PMNS mixing matrix. 
This allows us to accommodate sin^ 023 < 1/2 as well as sin^ 023 >1/2. In other cases 
the discussion of this permutation is already implicitly included in our analysis. In Case 
I and Case II we also confirm the estimate made for the solar mixing angle sin^ 0 i 2 that 
is bounded from below sin^ 0 i 2 ^ 0.276, while the lower bound in Case III is 1/3. The 
upper bound sin^ 0 i 2 < 0.345 found in Case IV-Pl and IV-P2 is obeyed as well, see table 
[T] and figure [TJ 

As can be read off from table j^for Case II and III, Ay^ = 10 “ ^min NO ~ 

showing that NO is better compatible with maximal atmospheric mixing sin^ 023 = 1 / 2 . 
This is simply due to the asymmetric 1 a errors of the atmospheric mixing angle for NO 
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and 10, see Q. Concerning the results for Case IV-P2 one might have naively expected 
that X mm 10 smaller than X mm NO’ value obtained for the atmospheric 

mixing angle at the best htting point ~ 0.255 is larger than 7r/4, i.e. sin^ 023 > 1/2- 
However, due to the large la error associated with sin^ 023 in the case of NO, for a 
value of the latter larger than (sin^ 023 )no ~ 0.452, namely +0.052 and a smaller one in 
the case of 10, for values of sin^ 023 smaller than the best fit value (sin^ ^ 23)10 = 0.579, 
namely —0.037, see Q, respectively, we find that NO hts slightly better in this case. 

The results found in table can be nicely visualized in the different sin^ 0jj planes, 
see figures and We plot the experimentally preferred 1 cr, 2 u and 3 a areas at two 
degrees of freedom in different colors (green, yellow and blue, respectively) for a neutrino 
mass spectrum with NO and 10 (left and right panels, respectively) in the sin^ dij planes 
using the data sets available in In order not to penalize NO we subtract = 0.97 
in this case so that the minimum value of Ay^ is zero for both mass orderings. The 
experimental best fit values of sin^ 9ij are presented by a cross in each plane. The shaded 
areas are to be compared with the black curves parameterized with 0 shown for Case I 
through IV-P2. On these curves we indicate the point (s) at which the y^ function 
for NO and 10, respectively, is minimized with a red dot. Furthermore, an arrow on 
each curve marks the direction of increasing 0 in the interval 0 < 0 < 0bf(,i)- We use 
the convention that the arrow always belongs to the closest label I, II, etc. of a certain 
case. Since the curve representing the results of Case I appears to be disconnected in 
the upper plots in figure due to the chosen scales of the axes we indicate different 
values of the parameter 0 with different symbols on the curve. In figure dashed red 
vertical lines show the lower and upper bounds on sin^ 0 i 2 that exist in the various cases. 


see (31), and comment below (57). Several curves (partly) overlap in the different 


planes: in figure and the upper two plots in figure clearly the curves belonging to 
Case II and III overlap, since both these cases predict maximal atmospheric mixing; 
furthermore, in the lower plots in figure the curves of Case I and II as well as of Case 
IV-Pl and Case IV-P2 lie on top of each other, since the relation between the solar and 
the reactor mixing angles is identical in these cases. Eventually, we note that in the 
plots showing sin^ 023 on the vertical axis in figures and the fact that the patterns 
of Case IV-Pl and IV-P2 are related by the exchange of the second and third rows in 
the PMNS mixing matrix is clearly visible, since their curves are symmetric with respect 
to the value sin^ 023 = 1/2. The curve belonging to Case I itself possesses this property 
showing that values of sin^ 023 smaller or larger than 1 /2 can be achieved for different 
choices of the parameter 0. For this reason the best fitting point 0bf is at a (very) 
different position for NO and for 10 in figure as well as in the upper plots of figure 
This does not occur in the other cases. 

We note that for Case I there exists a second (local) minimum of the y^ function for 
NO and 10 with y^; < 27. As can be guessed this second minimum for NO is obtained 


Those for sin^ 0i3-sin^ 9\2 and sin^ 0i3-sin^ 023 are explicitly given in [^, while the one for the third 
sin^ Oij plane can be simply constructed by summing up the values of the one-dimensional y^ projections 
for sin^ 012 and sin^ 023. This is justified, since these two mixing angles are to good approximation 
uncorrelated. We thank Thomas Schwetz-Mangold for help regarding this point. 
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1.8 2.0 2.2 2.4 2.6 


sin^0i3/lO ^ 




Figure 2: Results for the mixing angles sin^ ^ 23,12 with respect to sin^ 6*13 for Case I to III and 
Case IV-Pl and IV-P2. For conventions see caption of figure [7J The two curves for Case I in 
the upper two plots appear not to be connected due to the range plotted. Furthermore, note that 
the curves belonging to Case II and Case III overlap in the sin^ 0 i 3 -sin^ 023 plane. Similarly, the 
curves for Case I and II as well as for Case IV-Pl and IV-P2, respectively, lie on top of each 
other in the sin^ 0 i 3 -sin^ 0 i 2 plane. 


at ~ 2.967 for x^in ~ 10-42 and the mixing angles are very similar to those achieved 
at the global minimum of 10, sin^ 612 ~ 0.283, sin^ 0 i 3 ss 0.0217 and sin^ 6*23 ~ 0.592. In 
the same vein we find a second (local) minimum of the function for 10 at the best 
fitting point ~ 0.174 for x^in ~ 24.52. The values obtained for the mixing angles are 
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practically those of the global minimum of NO, i.e. sin^ 0i2 ~ 0.283, sin^ 0i3 ~ 0.0217 
and sin^ 023 ~ 0.408, compare table In all other cases the global minimum of the 
function that we mention in table is the only minimum with X^in ~ 

A closer look at table reveals that the solar mixing angle differs by up to 20% 
between Case I/Case II with sin^ 0i2 ~ 0.283 and Case Ill/Case IV-PI and IV-P2 where 
larger values of sin^ 0i2 are obtained. The experiment JUNO 29 will be able to reduce 
the error on the best fit value of sin^ 0i2 to ~ 0.7% at the 2 a level, thus allowing for a 
discrimination among Case I/Case II and Case Ill/Case IV-PI and IV-P2. According to 
the RENO-50 collaboration their planned experiment can achieve a similar reduction of 
the error [^. For the atmospheric mixing angle, no distinction is possible among Case 
II and III, since 023 is maximal in both cases. However, the predictions for Case I and 
Case IV-Pl/Case IV-P2 considerably differ: for NO (10) the atmospheric mixing angle 
is smaller (larger) in Case I than in Case IV-PI (Case IV-P2). This difference is large 
enough to be possibly distinguished in the experiment NOi^A [^. This experiment can 
also help in measuring the Dirac phase 5, especially, if running in both, the neutrino and 
the anti-neutrino, modes so that a discrimination between Case Il/Case III where <5 is 
maximal and Case I/Case IV-PI and IV-P2 with J = 0, vr might be possible. In contrast, 
the predictions for 0i3 are almost the same in Case I through Case IV-PI and IV-P2 
and, hence, it is unlikely that they can be distinguished at future neutrino facilities. 

Furthermore, we mention the outcome of our x^ analysis for the additional case for 
Ge = where we can choose as representative generators of the residual symmetries 
{Q,Z,X) = (T,S,Xo). We remind the reader that this case fails in giving a good 
fit to the solar mixing angle. We find two best fitting points i ~ 0.283(4) and 
0 bf2 ~ 2.858 for NO (10) whose sum approximately equals tt. Like in the cases above, 
this is expected from some symmetry of the formulae for the solar and the reactor mixing 
angles. The minimal x^ values for NO and 10 are at these two points X^in NO ~ 14.08 
and Xmin lO ~ 17-83, respectively. Thus, also here Ax^, the difference of Xmin 10 
X mm NO’ about 3.7, since the atmospheric mixing angle is fixed to be maximal by 
this pattern. The values obtained for the other mixing angles are sin^ 0i2 ~ 0.260 and 
sin^ 013 ~ 0.0216(7) and the Dirac phase 6 is maximal. In particular, we find for 0 = 0bf i, 
sin J = —I and Jcp ~ —0.0315(6) for NO (10) and for the other value 0 = 0bf,2) sinJ = 1 
and a positive sign also for Jcp- Like in the cases presented in table the two Majorana 
phases a and /3 are trivial. 

All patterns that can accommodate the experimental data well for a certain value 
of 0 predict trivial Majorana phases independently of 0. However, this is not a general 
result of the scenario with the flavor group and a CP symmetry. Instead there are 
also patterns that lead to three non-trivial CP phases which depend on the parameter 0 
and that are also non-trivial at the best fitting point(s) 0bf. For Ge being a Klein group, 
there exists a further category which can accommodate three non-trivial mixing angles. 


see discussion at the beginning of subsection 3.3 This category also predicts non-trivial 
and non-maximal Dirac and Majorana phases. A representative is {{Qi,Q 2 } ,Z,X) = 
{{S,T^ST^ST‘^} ,ST^ST‘^S, STXo). The third column of the PMNS mixing matrix 
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has then components that are the same as in (51) up to permutations. Thus, the reactor 
and the atmospheric mixing angles do not depend on 9 and are outside the experimentally 
preferred 3cr ranges. As a consequence, the values for NO and 10 are (much) larger 
than 100 and it appears in general difficult to achieve good agreement with the data, 
even if corrections to the leading order results presented here are expected to exist in an 
explicit model. 


4 Summary 

We have discussed a scenario for three Majorana neutrinos with the flavor group 
and a CP symmetry. These symmetries are broken to residual groups Ge = Z^, or 
Z2 X Z2 and to Gy = Z2 X CP the charged lepton and neutrino sectors, respectively. 
As a consequence, lepton mixing angles as well as CP phases can be predicted in terms 
of a single free parameter 9 that varies between 0 and vr. We have comprehensively 
studied all possible residual groups Ge, possible choices of Z2 being a subgroup of A 5 and 
CP symmetries that can be consistently combined. We have shown explicitly that the 
automorphisms corresponding to the CP transformations we consider are class-inverting 
and involutive and, furthermore, that these are the only automorphisms of A 5 with such 
properties. Performing a detailed analytical and numerical study, we found that only four 
mixing patterns exist (two for Ge = Z^ and one each for the other two choices, Ge = -^3 
and Ge = Z2X Z2) that can accommodate the experimental data on lepton mixing angles 
at the 3 a level or better for a particular value of the parameter 9. The reactor mixing 
angle is usually accommodated very well, while the solar one is bounded from below by 
three of the four patterns (sin^ 6*12 > 0.276 for Case I and Case II and sin^ 9 i2 > 1/3 for 
Case III, since this pattern has a TM column) and bounded from above, sin^ 9 i2 < 0.345, 
by the pattern called Case IV-Pl/Case IV-P2. Interestingly enough, two of the four 
patterns, Case II and Case III, predict a maximal Dirac phase 6 (together with maximal 
atmospheric mixing), while the other two patterns lead to sin<5 = 0 (and @23 is in 
general non-maximal). Common to all patterns is the prediction of trivial Majorana 
phases. Thus, in two of the four setups an accidental CP symmetry, common to the 
charged lepton and neutrino sectors, is present. Furthermore, we mention the existence 
of one further case for the choice Ge = Z^ that can be obtained by using as generators of 
the residual symmetries (Q, Z, X) = (T, S, Xq). This case fits the experimental data well 
to a certain extent. It mainly fails because of a too small value of the solar mixing angle, 
sin^ 9 i2 ~ 0.260. However, this could be easily reconciled with the experimental results 
in an explicit model with small corrections. The further predictions of this pattern are 
practically identical to those of Case II, i.e. the Dirac phase and the atmospheric mixing 
angle are maximal, the two Majorana phases are trivial and the reactor mixing angle can 
be accommodated well. As mentioned, all patterns that are preferred by experimental 
data turn out to lead to trivial Majorana phases. However, this is not a generic feature 
of our scenario with the flavor group A 5 and a CP symmetry, but occurs, because none 
of the patterns with non-trivial Majorana phases admits a reasonable fit to the data on 
lepton mixing angles. Interestingly enough, the absence of non-trivial Majorana phases 
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and the correlation between having a maximal Dirac phase and maximal atmospheric 
mixing or a trivial Dirac phase and non-maximal 623 for patterns that fit the data well 
have already been observed in analyses of scenarios with the flavor group Gf = S4 (or 
Gf = A4) and a CP symmetry j9, 
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Note added: At the very final stages of the completion of this work a paper 32 


dealing with the same topic appeared on the arXiv. Our results agree with those obtained 
However, the analysis in this work has been performed using results of a 


m 


32 


different global fit [^. In addition, the authors of display results for neutrinoless 
double beta decay and some model realization. 
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A Comments on automorphisms of A 5 


The CP transformations in ( |16[ ) all correspond to inner automorphisms, see ( |20[ ). We can 
check that these are class-inverting 11 and involutive. The latter feature is required, 
since X has to be a symmetric matrix in the flavor space, as explained below Q. Both 
properties are in particular given, if the twisted Frobenius-Schur indicator et(r) of an 
automorphism b is -|-1 in all irreps. r. The definition of ei(r) in the case at hand is 


60 


Xrig'g) 


(65) 


96A5 


with Xr{h) being the character of the element h and ''g the image of the element g under 
the automorphism i. An easy way to check 
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equality 




^Xr(e) = |{ff G As I'-£? = £/ ^}| 


that et(r) = -|-1 for all r is to verify the 

( 66 ) 


with Xr(e) being the character of the neutral element e in the irrep r, i.e. we sum 


over the dimensions of all irreps. of A 5 on the left-hand side of ( 66 ). For the trivial 


automorphism, see (19), the equality in ( 66 ) is obvious: the left-hand side is the sum 


of the dimensions of the irreps. r that is 16, while the right-hand side is the number of 
elements of the group A 5 that fulfill g = g~^. Clearly these are all Z2 generating elements 
together with the neutral element. Thus, these are 16 elements. For the automorphisms 


corresponding to choices of X other than Aq, see (16), we can do this computation using 
GAP ( 34 ] and also verify in these cases that the twisted Frobenius-Schur indicator for 
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all irreps r is +1. Thus, all CP transformations in (16) correspond to class-inverting, 


involutive automorphisms. Furthermore, we can compute the twisted Frobenius-Schur 
indicator for all other automorphisms i and in all representations r of and find that 
none of the other ones fulfills et(r) = ±1 for all r. So, none of the other automorphisms 
of is class-inverting and involutive. 

B Convention for lepton mixing parameters and CP in¬ 
variants 


We use the following convention for the PMNS mixing matrix 

Upmns = C/diag(l,e'“/2^e'(^/2+5)^ 

with U being defined, similar to the Cabibbo-Kobayashi-Maskawa mixing matrix 
0 


(67) 
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U = 


1 0 
0 cos 023 sin 023 
0 - sin 023 cos 023 


cos 013 0 sin 013 6 

0 10 

— sin 013 6®*^ 0 cos 01 
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cos 012 sin 012 0 
— sin 012 cos 012 0 
0 0 1 
( 68 ) 

The mixing angles Oij are taken to be in the interval between 0 and tt/2. The Dirac 
phase 6 as well as the two Majorana phases a and /3 can assume values between 0 and 


2 TT. The Dirac phase 5 can be extracted using the Jarlskog invariant Jcp 36 


JCP = ^^{UpMNS,llUpMNS,l^^PMNS,'!il^PMNS,'ii) 

= - sin 2012 sin 2023 sin 20i3 cos 0i3 sin 5 . 

8 

Similar invariants, called Ii and I2, can be defined for the Majorana phases 

Ii = Im i2(C7pjifiY5 ii)^) = sin^ 012 cos^ 012 cos^ 013 sina 

h = Im(t/|.;if ^5 i 3 (t/p^^ 5 ^id^) =sin^ 013 cos^ 012 cos^ 013 sin/3 


(69) 


(70) 

( 71 ) 


Notice that the Dirac phase 5 has a physical meaning only if all mixing angles are 
different from 0 and 7r/2. Analogously, the vanishing of the invariants Ii ^2 only implies 
sin a = 0, sin/3 = 0, if solutions with sin20i2 = 0, cos0i3 = 0 or sin20i3 = 0, cos0i2 = 0 
are discarded. Furthermore, notice that one of the Majorana phases becomes unphysical, 
if the lightest neutrino mass vanishes. 
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